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We extend the QCD Lagrangian to include Pauli-Villars (PV) gluons, quarks, and ghosts in such 
a way as to retain BRST invariance in an arbitrary covariant gauge. The extended Lagrangian 
can provide a starting point for nonperturbative calculations in QCD, particularly with light-front 
techniques, and the methods used to construct it may be useful for perturbative calculations in 
theories where dimensional regularization is not viable. The regularization is arranged by having 
all interaction terms in the Lagrangian be couplings between null fields, specific combinations of 
positive and negative-metric PV fields. The construction is done in steps, beginning with a gauge- 
invariant Lagrangian with massless PV gluons and degenerate-mass PV quarks. Auxiliary scalars 
are introduced to give mass to the PV gluons, and break the mass degeneracy of the PV quarks, 
following a method due to Stueckelberg. Gauge fixing terms for the gluon fields are a part of this 
construction. The ghost terms are then obtained and shown to provide the BRST invariance. A 
lack of dependence on the gauge parameter can be checked by the calculation of physical quantities 
for a range of values of the parameter. 


The principal goal of hadronic physics is the nonper¬ 
turbative solution of quantum chrornodynamics (QCD), 
in order to compute properties of hadrons as bound states 
of quarks and gluons. Lattice gauge theory [l] is one ap¬ 
proach for doing this; it has met with considerable suc¬ 
cess, but is limited by its reliance on a Euclidean frame¬ 
work and by other difficulties, such as fermion doubling 
and a large pion mass. Dyson-Schwinger methods 0 are 
also useful but also Euclidean. The truncated conformal 
space approach {§[ shows promise but is in its infancy. 
An alternative, which retains a Minkowski framework 
and the intuitive notion of wave functions, is the light- 
front Hamiltonian approach Si- Light-front methods 
have been successfully applied to many gauge field the¬ 
ories, including quantum electrodynamics [6|-l8|, QCD 
in two dimensions Q, and supersymmetric Yang-Mills 
theories [10]; however, with respect to four-dimensional 
QCD, this approach has always lacked a consistent reg¬ 
ularization that can be used nonperturbatively. 

The standard regularization used for non-Abelian the¬ 
ories such as QCD is dimensional regularization 11]. Un¬ 
fortunately, this method is inherently perturbative, rely¬ 
ing as it does on the analysis of the singularity structure 
of Feynman diagrams and the modification of the asso¬ 
ciated integrations. All of this is buried deep within a 
nonperturbative calculation, where such separations and 
modifications are difficult if not impossible to perform. 

A regularization that can be introduced at the La¬ 
grangian or Hamiltonian level is much to be preferred, 
and there has been just such an approach for Abelian 
theories, developed in the context of light-front Hamil¬ 
tonian methods dm, 00 Massive Pauli Villars (PV) 


fields are introduced to the Lagrangian with couplings 
arranged to provide the cancellations necessary for reg¬ 
ularization via subtractions of negative-metric PV con¬ 
tributions, and the Hamiltonian derived in some fixed 
gauge. 

The cancellations are obtained by constructing all in¬ 
teractions from combinations of positive and negative- 
metric fields that are null in the following sense. Let the 
(light-front) mode expansions for the gluon and quark 
fields be written as 
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with dk = dk + d 2 k± the light-front volume element [||, 
a set of polarization vectors, a\ kx the gluon creation 
operator with color index a, PV index k , and polarization 
A, and b' ais the quark creation operator with color index 
a, PV index i, and spin s. The quark flavor index is sup¬ 
pressed, as it will be throughout. The physical fields are 
denoted by a PV index of zero. The (anti)commutation 
relations for the creation operators are 


Wakx(q),al lx ,(i j')] =r k e x 6 ab 5 k i6\y6(q-(/), (3) 
{bais(q),blj S 'W)} = SiSabSijS SS 'S(q- q). (4) 


1 There has also been work on construction of a light-front regu¬ 


larization that is perturbatively equivalent to covariant regular¬ 
ization M- 
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Here r/. = ±1 and Sj = ±1 are metric signatures for the 
gluon and quark fields with PV index k and i, respec¬ 
tively. In addition, the gluon field has a polarization de¬ 
pendent metric signature specified by e A = (—1,1,1,1). 
Null fields are then constructed as 

= (5) 

k i 

where the and are coupling coefficients constrained 
to satisfy 

E r ^ = 0 ’ E s ‘A 2 = o. (6) 

k i 

Obviously, at least one PV gluon and one PV quark must 
have a negative metric. Typically one would take rk = 
(-l) fc and Si = (-l)L 

Interaction terms, such as gij>'y IJ, T a A aki iib, are then built 
from null combinations, with rg = sq = £o = A) = 1 for 
the physical field. Here T a is a matrix representing the 
Lie algebra, and repeated color indices are summed. Any 
loop must then contain the sum Ylik s iPiPi r ktikt; k , with 
the primed coefficients coming from the second vertex, 
which may or may not be the same as the first. The met¬ 
ric signatures come from contraction of the creation and 
annihilation operators associated with the internal lines 
of the loop. If the two vertices are the same, the null 
constraints immediately provide two subtractions; if not 
the same, the different null combinations from the differ¬ 
ent vertices must be mutually null, that is ]Th Sj/?j/3' = 0 
and = 0; to again provide two subtractions. 

This approach has been used to study QED in an arbi¬ 
trary covariant gauge, including study of the dependence 
of physical quantities on the gauge parameter {§]. The 
couplings induce currents that change the PV index of 
the held, including mixing with the physical held, which 
breaks gauge invariance; however, in the limit that the 
PV helds are removed, gauge invariance is restored. 

The extension to non-Abelian theories is, of course, 
problematic. Unlike Abelian theories, where a massive 
vector boson is known to be renormalizable fl5j | , the proof 
of renormalizability for Yang-Mills theories is best ap¬ 
proached in terms of BRST invariance [lj|, which would 
appear to require an underlying gauge invariance, which 
in turn requires massless vector bosons. The well-known 
exception is the use of spontaneous symmetry breaking to 
provide mass to gauge bosons without destroying gauge 
invariance at the Lagrangian level, which is a hint as to 
how one might proceed. Another difficulty is that, even 
for massless PV gluons and mass-degenerate PV quarks, 
the null couplings that mix fields do break gauge invari¬ 
ance. 

We have resolved these difficulties. The giving of 
mass to the PV gluons and the lifting of the PV-quark 
masses is accomplished by a non-Abelian generalization 
of Stueckelberg’s mechanism , which associates 


an auxiliary real scalar with each gluon field, simulta¬ 
neously giving mass to the gluon and the scalar while 
also fixing the gauge. The breaking of gauge invariance 
by the field-mixing interactions is eliminated by general¬ 
izing the definition of the gauge transformation to also 
include field mixing; the original gauge transformation is 
recovered in the limit that the PV fields are removed, by 
taking their masses to infinity. The BRST invariance is 
recovered for finite PV masses by inclusion of the appro¬ 
priate Faddeev- Popov [l9j] ghost terms required by the 
gauge-fixing terms. 

In case this sounds too good to be true, we must point 
out that there are drawbacks to our approach. One is the 
presence of nonlocal interactions for the ghost fields; this 
is made necessary by the mechanism used to generate PV 
gluon masses. The other drawback is a proliferation of 
PV fields, which will be a large computational burden. 
There are several constraints that the PV couplings must 
satisfy, and each constraint requires the addition of one 
or two PV fields. The formulation given here can require 
as many as four PV gluons, three PV quarks, five PV 
adjoint real scalars, and four PV ghosts and anti-ghosts. 
This is for each color and flavor, and therefore is a very 
large number. 

Oddly enough, the Higgs mechanism is not used. We 
found that the requirement of null couplings, to provide 
the regularization, causes two difficulties. One is in the 
Higgs sector itself, where the null (fA interaction does not 
lead to a well-defined minimum to drive the breaking of 
the symmetry. The other is that gauge-invariant null cou¬ 
plings of the PV Higgs to the gluons requires that any 
symmetry breaking leave as massless only a null com¬ 
bination of gluon fields, rather than produce a massless 
physical gluon and massive PV gluons. 

In what follows, we describe the construction in stages. 
The first is a Lagrangian for massless gluons and mass- 
degenerate quarks that is invariant with respect to a gen¬ 
eralized gauge transformation. Next, we introduce the 
auxiliary scalars and add terms that give mass to the PV 
gluons and fix the gauge. This is followed by specifica¬ 
tion of a term that lifts the mass degeneracy of the PV 
quarks. Finally, we determine the Faddeev-Popov ghost 
term and state the BRST transformations of the fields. 
Our construction is kept slightly more general than QCD, 
in that we treat SU(N) Yang-Mills theory with funda¬ 
mental matter, with QCD being the N = 3 case. 

The gauge-invariant Lagrangian for massless PV glu¬ 
ons and mass-degenerate PV quarks is 

£ = -J E r k F ak F aknv + E - 1TI )^i 

k i 

+9 E PiPj^il^TgAak^j. (7) 

ijk 

The quark spinor field t[>i is a column vector with respect 
to color. The quarks (of a single flavor) are all of mass 
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m. The gluon field tensor F^ k is computed from the field 
as 

Kk = ^Kk - d V A: k - r^gfabc E (8) 

Im 

The structure constants f a bc specify the commutation re¬ 
lation [T a ,T b ] = ifabcT c . 

The gauge transformations of the fields are 

^ak —* A^ k + <9 M A a k + r k ^ k gfabcA b A^ (9) 

A —> A + igsi/3iT a A a ip, (10) 

with A a = ^2 k £ k A a k- The null combination ifj is then 
gauge invariant, and the null combination A% is Abelian 

AZ-+AZ + &* A a , (11) 


as is the associated field tensor 


k 


( 12 ) 


The Lagrangian C is gauge-invariant with respect to 
these transformations. 

When expressed in terms of the null combinations, the 
Lagrangian reduces to 

c = E - d v K k ? + gfabcd^A-A^ 

k 

+ E - m)A + gi’l^TaAa^. (13) 


All interactions are between null combinations, and 
the free part of the Lagrangian corresponds to terms 
for massless gluon fields with metric signatures r k and 
degenerate-mass quark fields with signatures .s t . The 
four-gluon interaction that normally appears in the QCD 
Lagrangian has disappeared, as a result of cancellations 
between physical and PV gluons. It is restored between 
physical gluons in the infinite-mass limit for PV gluons 
as a contraction of two three-gluon vertices, with the con¬ 
traction being a PV gluon. This provides a mechanism, 
natural in the context of the present Lagrangian, for the 
often-used trick of introducing an auxiliary field to reduce 
the four-gluon interaction to two three-gluon interactions 
for the convenience of color factors in perturbation the¬ 
ory [Io|]. 

To give mass to the PV gluons, we use a nomAbelian 
extension of the Stueckelberg mechanism 17, 18]. Real 
adjoint scalar fields cj) ak are introduced with a gauge 
transformation of 

( fiak ^ (frak “b [A k A(lk T nk^k^kgfabc I dx gA b (x )A(((x ). 

04 ) 

The line integral is present in order that the gauge trans¬ 
formation of the derivative be 


d^Ak —> d^fiak + iikd^Aak + nkrAkgfabcA b A^. (15) 


This then makes the combination g- k A^ k — d^Ak gauge 
invariant, which is used in the construction of the first 
term of an additional piece for the Lagrangian, written 
as 

A = \ E Tk ~ d'Vafcf (16) 

k 

E rfc ( 9 ^afc + ^Ak'j • 


The second term is the gauge fixing term. However, when 
the two terms are combined, the cross terms sum to a 
total divergence which can be neglected, leaving 


A 


\ E rk ^k ( A akf - | E rk ( 8 K A akf 



iAAkf 


/4 



(17) 


This is immediately seen to provide a mass g k and stan¬ 
dard gauge-fixing term for each gluon field, with gauge 
parameter and a free Lagrangian for real scalars with 
masses g k /y/~C, and metric signatures r k . For the physical 
gluon, the mass no is to be taken to zero. For pure Yang- 
Mills theory, this can be done explicitly, by exclusion of 
the scalar fields <f> a o and replacement of the k = 0 terms 
in C g by a simple gauge-fixing term for the physical glu¬ 
ons: — | (i9 m A(] 0 ) . When quarks are included, the k = 0 
scalar is needed, and the mass no must be taken to zero 
as a limit; however, as will be seen, the scalar <j > a o does 
not couple to the physical fields for any value of no- 
The mass degeneracy of the PV quarks can be re¬ 
moved by coupling the quarks to a null combination of 
the scalars <j> a k, 


A = E£ fc ^7“ fe ’ ( 18 ) 

fc 

made null by the additional constraint 

e^ 4=°. (i9) 

k 

Here /ipy = max*, nk is the mass scale of the PV gluons. 
The gauge transformation of the combination is Abelian 

4*a -> 4>a + nPvAa (20) 

and provides the necessary piece to make gauge-invariant 
the Lagrangian term 

£q = - E Simi{A + ig-^-A$T a )(A - ig^-4> a T a ip), 

i npv gpv 

( 2 !) 

which is to replace the quark-mass term — 

in m- Thus, each PV quark can have a different mass 

rrii, with mo the mass of the physical quark. 
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This new term can be written as 

TOpy 


C q = -^2 - ig 


Mpv l 


1pT a (p a 1p - 1pT a (p a 1p 


~ S i m f l 1pT a (p a Tb(pb1p, 

i ^PV 


( 22 ) 


where we define a PV-quark mass scale mpy = max^ rrii 
and a new null combination of quark fields 


mpy 


(23) 


For this to be null, there is the additional constraint 

= 0, (24) 


number of PV gluons; the number of PV scalars is one 
more, except for pure Yang-Mills theory, where (p a o is 
not needed. Because fig -C /ipy, there could be difficulty 
with significant figures in numerical solutions of these 
constraints. However, calculations in QED [§} have indi¬ 
cated that the no //xpy ratio need not be extremely small, 
so that this potential difficulty may be mitigated. 

This completes the construction. The full expression 
for the PV-regulated QCD Lagrangian can be obtained 
by adding C g and £fp, found in (fT71) and (l26l) . to £ in 
m, and replacing the quark mass term in £ with £ g , 
found in (l22l) . 

The associated BRST transformations are determined 
from the gauge transformations with the replacement 
A a k —> eCa/c, where e is a real Grassmann constant for 
which e 2 = 0. The transformations are 


and for ip and ip to be mutually null, we must have 


^SirriiPt = 0. (25) 


This last constraint has the advantage of making the 
last term in (OS simply zero. Also, the second and 
third terms combined are proportional to ■ /3i/3j(rrij — 

rrii)ipiT a (p a ipj , which is zero for i = j and guarantees 
that none of the scalars are directly coupled to the phys¬ 
ical quarks, as claimed above. There are now three con¬ 
straints on the PV quark masses and coupling coeffi¬ 
cients, which requires three PV quarks if all masses m; 
are to be chosen independently. 

The last step of the construction is the Faddeev-Popov 
ghost term, which can be done by the usual methods of 
path-integral quantization [2l| . We obtain, for ghosts c a k 
and anti-ghosts c a k, 


„ 2 

£fP — ^ ^ Xkd^Cakd^ Cak ^ ^ Pfc ~ ^-ak^ak (26) 

k k ’ 


~^~g fabc 


d^CaCbA^ - ^r-c a J dx^c b (x')A%(x') 


with the null combinations defined as 


Mfe 


— / — y ^kX-ak: C a — / o ^ak- 

^ ^ ^ Mp V 


k k k 

For these to be (mutually) null, we require 


(27) 


r k^Uk = 0, r ^Uk = 0- (28) 

k k 

There are now four constraints on the PV-gluon masses 
Hk and coupling coefficients the two above plus J6]) 
and (fTiil) . If all the PV-gluon masses are to be chosen 
independently, this would require four PV gluons; for 
constrained mass values, the number of PV gluons can 
be less. The number of PV ghosts is the same as the 


SA ak = ed ^ c ak + tVk^kgfabcCbA (29) 
Hi = iegSiPiT a c a Hi = -iegSiPi^TaCa, (30) 
fifiak ^l^k^ak (31) 


Aerkik^kgfabc 



c a k + er k ^ k fi k gf abc c b A^ (32) 

Sc ak = — C e (dn A ak + y<t>a.kj , (33) 

dc a k — -er k £kgfabcCbC c . ( 34 ) 


For the various null combinations we then find 


6A% = ed^Ca, S(p a = e/rpy c a , Sc a = 0, (35) 

dip = 0, Sip = 0. (36) 


All of the gauge-invariant pieces of the Lagrangian are 
then automatically BRST invariant, and we need only 
check the sum of the second term of £ g and £fp , in (fT6l) 
and (l26l) respectively. A short calculation shows that this 
sum, and therefore the entire Lagrangian, is also BRST 
invariant. 

However, what we have done is only a formal construc¬ 
tion. A next step is to construct the corresponding light- 
front Fock-space Hamiltonian and then to solve for its 
eigenstates. This will be facilitated by another aspect of 
null couplings, in that the projections 0 ip- = ^7°7 ip 
and ip- = ^7°7 ~ip of the Dirac field are gauge invari¬ 
ant and satisfy constraint equations, projected from the 
Dirac equation, that are independent of the gluon fields. 
This allows for explicit solution of the constraint with¬ 
out invocation of light-cone gauge [Q]. There will also be 
cancellation of instantaneous fermion interactions, just 
as seen in Yukawa theory [l2j and QED This should 
be a promising avenue of attack on QCD. 

At the stage of solving the Fock-space eigenvalue prob¬ 
lems, there is a potential source of gauge dependence, 
in that most methods employ Fock-space truncation to 
produce a finite computational problem. A method that 
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avoids such a truncation, the light-front coupled-cluster 
method [22], has been developed and should be applica¬ 
ble to QCD. Rather than truncate Fock space, a finite 
computational problem is obtained by truncation of the 
way in which higher Fock-state wave functions are related 
to the lower ones. 


[1] C. Gattringer and C.B. Lang, Quantum Chromodynamics 
on the Lattice, (Springer, Berlin, 2010); H. Rothe, Lattice 
Gauge Theories: An Introduction, 4e, (World Scientific, 
Singapore, 2012). 

[2] I. C. Cloet and C. D. Roberts, Prog. Part. Nucl. Phys. 
77, 1 (2014). 

[3] M. Hogervorst, S. Rychkov and B. C. van Rees, Phys. 
Rev. D 91, 025005 (2015). 

[4] S.J. Brodsky, H.-C. Pauli, and S.S. Pinsky, Phys. Rep. 
301, 299 (1998); J. Carbonell, B. Desplanques, V.A. Kar¬ 
manov, and J.F. Mathiot, Phys. Rep. 300, 215 (1998); 

G. A. Miller, Prog. Part. Nucl. Pliys. 45, 83 (2000); M. 
Burkardt, Adv. Nucl. Phys. 23, 1 (2002); B. L. G. Bakker 
et al., Nucl. Phys. Proc. Suppl. 251-252, 165 (2014). 

[5] S. D. Glazek, Phys. Rev. D 85, 125018 (2012); 

M. Gomez-Rocha and S. D. Glazek, arXiv: 1505.06688 
[hep-ph], 

[6] S.J. Brodsky, V.A. Franke, J.R. Hiller, G. McCartor, S.A. 
Paston, and E.V. Prokhvatilov, Nucl. Phys. B 703, 333 
(2004); S.S. Chabysheva and J.R. Hiller, Phys. Rev. D 
81, 074030 (2010); Phys. Rev. D 82, 034004 (2010). 

[7] J.P. Vary et al, Phys. Rev. C 81, 035205 (2010); 

H. Honkanen, P. Maris, J.P. Vary, and S.J. Brodsky, 
Phys. Rev. Lett. 106, 061603 (2011); X. Zhao, H. Honka¬ 
nen, P. Maris, J.P. Vary and S.J. Brodsky, Phys. Lett. B 
737, 65 (2014). 

[8] S.S. Chabysheva and J.R. Hiller, Phys. Rev. D 84, 
034001 (2011). 

[9] K. Hornbostel, S. J. Brodsky and H. C. Pauli, Phys. Rev. 


D 41, 3814 (1990). 

[10] Y. Matsumura, N. Sakai, and T. Sakai, Phys. Rev. D 52, 
2446 (1995); O. Lunin and S. Pinsky, AIP Conf. Proc. 
494, 140 (1999); J.R. Hiller, S.S. Pinsky, Y. Proestos, 
N. Salwen, and U. Trittmann, Phys. Rev. D 76, 045008 
(2007) and references therein. 

[11] G. ’t Hooft and M. Veltman, Nucl. Pliys. B 44, 189 
(1972); C.G. Bollini and J.J. Giambiagi, Phys. Lett. B 
40, 566 (1972); J.F. Ashmore, Nuovo Cim. Lett. 4, 289 
(1972). 

[12] S.J. Brodsky, J.R. Hiller, and G. McCartor, Ann. Phys. 
305, 266 (2003); 321, 1240 (2006). 

[13] V.A. Karmanov, J.-F. Mathiot, and A.V. Smirnov, Phys. 
Rev. D 77, 085028 (2008); Phys. Rev. D 82, 056010 
( 2010 ). 

[14] S.A. Paston and V.A. Franke, Theor. Math. Phys. 112, 
1117 (1997) [Teor. Mat. Fiz. 112, 399 (1997)]; S.A. Pas¬ 
ton, V.A. Franke, and E.V. Prokhvatilov, Theor. Math. 
Phys. 120, 1164 (1999) [Teor. Mat. Fiz. 120, 417 (1999)]; 
M. Y. Malyshev, S. A. Paston, E. V. Prokhvatilov, 
R. A. Zubov and V. A. Franke, arXiv: 1504.07951 [hep- 
th]; arXiv: 1505.00272 [hep-th], 

[15] P.T. Matthews, Phys. Rev. 76, 1254 (1949); A. Salam, 
Nucl. Phys. 18, 681 (1960; S. Kamefuchi, Nucl. Phys. 18, 
691 (1960). 

[16] C. Becchi, A. Rouet, and R. Stora, Phys. Lett. B 52, 
344 (1974); Commun. Math. Phys. 42, 127 (1975); Ann. 
Phys. 98, 287 (1976); I.V. Tyutin, Lebedev preprint 
FIAN No. 39, 1975, unpublished. 

[17] E. C. G. Stueckelberg, Helv. Phys. Acta 11, 225-244 
(1938); 11, 299-312 (1938); 11, 312-328 (1938). 

[18] R. Marnelius and I. S. Sogami, Int. J. Mod. Phys. A 13, 
3101 (1998). 

[19] L.D. Faddeev and V.N. Popov, Phys. Lett. B 25, 29 
(1967). 

[20] A. Grozin, Lectures on QED and QCD, (World Scientific, 
Singapore, 2007). 

[21] D. Bailin and A. Love, Introduction to Gauge Field The¬ 
ory, revised ed., (Institute of Physics, Bristol, 1993). 

[22] S.S. Chabysheva and J.R. Hiller, Phys. Lett. B 711, 417 

( 2012 ). 




